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We present preliminary results of a new lattice computation of hadronic matrix elements of baryon number 
violating operators which appear in the low-energy effective Lagrangian of (SUSY-)Grand Unified Theories. The 
contribution of irrelevant form factor which has caused an underestimate of the matrix elements in previous 
studies is subtracted in this calculation. Our results are 2~4 times larger than the most conservative values often 
employed in phenomenological analyses of nucleon decay with specific GUT models. 



1. Introduction 

Nucleon decay is an important consequence of 
Grand Unified Theories (GUTs), which can give 
strong constraints on such theories. Unfortu- 
nately, theoretical estimations of the lifetime and 
branching ratios suffer from large uncertainties. 
One of the sources of the uncertainties is the value 
of hadronic matrix elements from the nucleon to 
pseudo scalar meson through baryon number vi- 
olating operators which appear in the low-energy 
effective Lagrangian of GUTs. 

Conventionally these matrix elements are esti- 
mated by rewriting them in terms of the three- 
quark annihilation amplitude (O|0|p) using the 
soft-pion theorem. There have been several lat- 
tice calculations for (O|0|p) [Q-^. They have not 
yielded a definitive result, however. Moreover, 
the validity of the soft-pion theorem for the nu- 
cleon decay process has not been clarified. 

A pioneering lattice work to calculate (7r°|0|p) 
directly was made in [g[. We also made an ef- 
fort to advance it in several fronts W . After this 
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study, we have found that previous lattice cal- 
culations for the matrix elements including ours 
had a subtle point which caused an underestimate 
of the matrix elements. The problem is that the 
previous results contain the contribution of an ir- 
relevant form factor, which vanishes after using 
the equation of motion of the out-going lepton 
and setting its mass to zero. This contribution 
can not be subtracted using our previous data, 
essentially because we did not measure all neces- 
sary components of the three-point function. 

In this article we report results of a new simu- 
lation in which we measure all the necessary com- 
ponents of the three-point correlation function to 
disentangle the relevant form factor from the ir- 
relevant one. 

As in our previous calculation, we evaluate ma- 
trix elements of all dimension six baryon num- 
ber violating operators classified according to 
SU{3) X SU{2) X [/(I) 11], so as to cover vari- 
ous GUT models and decay processes. Also the 
calculations are made at physical quark mass and 
physical momentum in order to take account of 
the N ^ K process which is expected to be the 



dominant mode for SUSY-GUTs. 

2. Calculational method 

We wish to calculate the matrix element of a 
class of three quark operators between the nu- 
cleon N and a pseudoscalar meson PS. Using 
Lorentz and parity invariance, one finds that the 
matrix elements are described by two form fac- 
tors, 



{PS\{q,CPn/Lq2)PRjLq3\N) 



Aiq' 



)Pr/lN 



B{q^)P^ 



R/L 



(1) 



where qi represents a quark field, g is a four mo- 
mentum of the out-going Icpton and N and Pr/l 
stands for the nucleon spinor and chiral projec- 
tion operator, respectively. 

Multiplying by the lepton spinor and using the 
equation of motion, one can see that the second 
term is proportional to lepton mass. We call A 
and B as relevant and irrelevant form factor, re- 
spectively, as the latter does not contribute when 
the lepton mass is neglected. 

To extract the form factors, we form a ratio, 

{Jps{ty)O^{t^)JN,-^'{0)) 



{Jps{ty)jls{tx)){JN{tx)JN{0)) 



where 7 and 7' are spinor indices, and Zps (Zn) 
is the residue of the two-point function of PS 
meson (nucleon). This ratio has the following 
asymptotic form, 



A + q"B 
Bq-a 



(3) 



in a 2 X 2 block notation. 

It is important that the upper component is a 
linear combination of the two form factors and the 
lower component is proportional to the irrelevant 
term. In order to disentangle the relevant form 
factor from the irrelevant one, both components 
are needed, while, in the previous studies, only 
the upper component was measured. It should 
be noted that we can not follow this procedure 
for the case q=0. 

To compare our results with the predictions of 
the tree-level chiral Lagrangian |g|,^ , we also cal- 
culate the three-quark annihilation amplitude, 



{0\{uCdR)uL\p) = aNL, 



(4) 



{0\{uCdL)uL\p) = PNl, (5) 

which are obtained from the two-point function. 

3. Numerical simulation 

Our calculation is carried out in quenched QCD 
at f3— 6.0 with the Wilson quark action on a 28^ x 
48 X 80 lattice. We analyze 100 configurations 
at the hopping parameter ii'=0. 15620, 0.15568, 
0.15516, 0.15464. The lattice scale fixed by 
mp=770 MeV in the chiral limit (^^=0. 15714(1)) 
is a^^=2.30(4) GeV, and the point for the strange 
quark estimated from ttik /mp=0.6AA is given by 
7^:^=0.15488(7). 

We fix the nucleon source at t = 0, PS meson 
sink at i = 29 and move the operator between 
them. Matrix elements are evaluated at four spa- 
tial momenta ga=(0,0,0), (7r/14,0,0), (0,7r/14,0), 
(0,0,7r/24), injected in the PS meson sink. As 
mentioned in the previous section, it is not pos- 
sible to disentangle the relevant form factor from 
the irrelevant one when ^'=0. Therefore, we use 
two data with two finite momenta to interpolate 
to the physical momentum. 

The form factor obtained by fitting the ratio 
of three-point function and two-point functions 
depends on the quark mass and the momentum 
squared. We distinguish u-d and s quark masses; 
the former {rriud) is taken to the chiral limit, 
and the latter (ms) interpolated to the physical s 
quark mass. We also interpolate the momentum 
square to the physical lepton mass. To do this, 
we fit the data with the following form, 

C1+C2 ■ q'^ +C3 ■ {q^Y + C4 • mud + C5 • m^. (6) 

Matrix elements are renormalized, with mixing 
included, by tadpole-improved one-loop renor- 
malization factors to the MS scheme calcu- 
lated at the scale /i = 1/a. The quoted errors are 
only statistical, estimated by the single elimina- 
tion jackknife procedure. 

4. Results 

We present a typical q^ dependence of the rele- 
vant form factor in Fig. ||. We also plot a lin- 
ear combination of the relevant and irrelevant 
form factors A + q^B and the prediction from 
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Figure 1. g^ dependence of the relevant form 
factor |ylTr(g^)| (circles) and a linear combina- 
tion of the relevant and irrelevant form factors 
l^niq"^) + q^ B.^{q'^)\ (up-triangles), obtained from 
the matrix element {TT'^\(uiiCdfj)uL\p) ■ The pre- 
diction from chiral Lagrangian is also shown 
(down-triangle) . 

tree-level chiral Lagrangian using |a|=0. 015(1) 
GeV'^ obtained in our simulation (we also obtain 
1/31=0.014(1) GeV^). We see that the value of 
the relevant form factor is close to that from the 
chiral Lagrangian. As expected from analytical 
considerations with the chiral Lagrangian, there 
is a cancellation between the relevant and irrel- 
evant form factors. We consider that this led 
to the inconsistency between the direct compu- 
tation through three-point functions and the in- 
direct computation using the tree-level chiral La- 
grangian reported in previous lattice studies MM ■ 
In Fig. 0, we plot our results for the relevant 
form factors at q^~0. To compare our results 
with the most conservative estimation of the ma- 
trix elements, often employed in phenomenologi- 
cal analyses of nucleon decay with specific GUT 
models, we also plot the prediction from tree-level 
chiral Lagrangian with a choice of the parameters 
|a| = |/3|=0.003 GeV^. Our results are about 2~4 
times larger than the most conservative values, 
implying a stronger constraint on the parameter 
space of GUT models. 
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Figure 2. Relevant form factors at physical mo- 
mentum g^ = from present calculation (filled 
circles) compared with the most conservative pre- 
dictions of tree-level chiral Lagrangian (aster- 
isks). 
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